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Recently there has been some controversy in the literature concerning the existence of a fermion 
sign problem in the jV — (2,2) supersymmetric Yang-Mills (SYM) theories on the lattice. In 
this work, we address this issue by conducting Monte Carlo simulations not only for JV — (2,2) 
but also for JV — (8,8) SYM in two dimensions for the U(N) theories with N = 2, using the 
new ideas derived from topological twisting followed by geometric discretization. Our results 
from simulations provide the evidence that these theories do not suffer from a sign problem as the 
continuum limit is approached. These results thus boost confidence that these new lattice formu- 
lations can be used successfully to explore the nonperturbative aspects of the four-dimensional 
,yV =4 SYM theory. 
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1. Introduction 

Supersymmetric Yang-Mills (SYM) theories form an interesting class of field theories from 
a variety of perspectives. They play a role as toy models for understanding the properties of non- 
abelian gauge theories including QCD. They could be strong candidates for describing the physics 
beyond the Standard Model. Some of these gauge theories also play an important role in the 
AdS/CFT correspondence, and thus are connected to a class of gravitational theories. Many inter- 
esting features are exhibited by these theories, for example, dynamical supersymmetry breaking, 
are inherently non-perturbative in nature. To study such non-perturbative features we are motivated 
to regularize them on a space-time lattice. 

Unfortunately, conventional regularization methods for supersymmetric field theories on the 
lattice have been proven difficult or even impossible in most cases. This difficulty can be traced 
back to the fact that the supersymmetry algebra, which is an extension of the usual Poincare algebra, 
is broken completely by naive discretization on a space-time lattice. Thus, one has to search for 
nontrivial methods to discretize these theories which respect at least part of the supersymmetry 
algebra. It turns out that Nature has provided us a set of elegant tools, topological twisting and 
orbifolding, to construct at least certain classes of supersymmetric theories on the lattice while 
preserving a subset of the continuum supersymmetry algebra. The reviews [|]] and references 
therein provide details of these techniques, which, we will be using in this work. There also exist 
other complementary approaches in the literature and they can be found in p]. 

Interestingly, several continuum SYM theories, including the well known JV = 4 SYM in 
four dimensions, can be implemented on a space-time lattice via geometric discretization of the 
corresponding topologically twisted forms of these theories. On a flat Euclidean space-time, these 
classes of continuum SYM theories and their topologically twisted cousins are related just through 
a change of field variables. 

Lattice theories constructed this way using these techniques have no doubling problem, respect 
gauge-invariance, preserve a subset of the original supersymmetries and target the usual continuum 
theories in the naive continuum limit. However, only certain classes of lattice SYM theories are 
possible with this construction scheme: the requirement is that the target SYM theories must pos- 
sess a sufficient number of extended supersymmetries. To be more precise, the number of super- 
charges must be an integer multiple of 2 D where D is the space-time dimension. This includes the 
Jf = (2, 2) SYM theory in two dimensions and JV = 4 SYM in four dimensions. In this work, we 
study these two theories, reduced in the case of the jV = 4 model to two dimensions, yielding the 
^ = (8,8) SYM theory. 

Even with the existence of a supersymmetric lattice construction for a given SYM theory, one 
might encounter another difficulty that would prevent one from extracting sensible results from lat- 
tice simulations. This difficulty is known as the fermionic sign problem. Consider a generic lattice 
field theory with a set of bosonic and fermionic degrees of freedom, <p and y respectively. Then 
the partition function of the theory can be written as Z = f{D<p][D\j/] exp(— Sb[<P] — y T M[<p]\i/). 
After integrating out the fermionic degrees of freedom, we have Z = f[D<p] Pf(M) exp(— Sb[0]). 
The matrix M corresponds to the fermion operator and Pf(M) the corresponding Pfaffian. For a 
2n x 2n matrix M, the Pfaffian is explicitly given as Pf(M) 2 = DetM. In the supersymmetric lat- 
tice constructions which we will consider in this work the matrix M at non zero lattice spacing 
is a complex operator, and one might worry that the resulting Pfaffian could exhibit a fluctuating 
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phase depending on the background boson fields . Since the Monte Carlo simulations must be 
performed with a positive definite measure, the only way to incorporate this fluctuating phase is 
through a reweighting procedure, which folds this phase in with the observables of the theory. 
Thus the expectation values of the observables derived from such simulations can be contaminated 
by drastic statistical errors, which could overwhelm the values of the quantities we are trying to 
measure. 

In the constructions of supersymmetric lattice gauge theories, there has been an ongoing de- 
bate on the existence of a sign problem in the two-dimensional j¥ = (2,2) lattice SYM [[| |], BJ]. 
In it was shown that there is a potential sign problem in the two-dimensional jV = (2,2) 
lattice SYM. Furthermore, in [Q] numerical evidence was presented of a sign problem in a phase 
quenched dynamical simulation of the theory at non-zero lattice spacing. More recently Hanada et 
al. [j|] have argued that there is no sign problem for this lattice theory as the continuum limit is 
approached. However, the models studied by these various groups differed in detail; Catterall et 
al. studied an SU (2) model obtained by truncating the supersymmetric U (2) theory and utilized 
bosonic link fields valued in the group SL(2,C), while Hanada et al. used a U (2) model where the 
complexified bosonic variables take their values in the algebra of U (2) together with the inclusion 
of supplementary mass terms to control the fluctuations of the scalar fields in the theory. 

In this work, we present results from simulations of the two dimensional j¥ = (2,2) U(2) 
SYM theory (which we will refer to from now on as the £2 = 4 theory, with £2 the number of 
supercharges) and the maximally supersymmetric jV = (8, 8) U (N) SYM theory (we refer to this 
theory as the £2 = 16 theory) using algebra based parameterizations of the bosonic link variables. 
We measure the average phase fluctuations as we approach the continuum limit. In all cases we 
find evidence that the phase fluctuations disappear in the continuum limit indicating the absence of 
the sign problem in these lattice theories as the continuum limit is approached. 

2. Supersymmetric Yang-Mills theories on the lattice 

As mentioned before we will be studying the lattice SYM theories constructed using the 
twisted approach followed by geometric discretization. At this point we also stress that the lattice 
SYM theories constructed using orbifold methods are equivalent to twisted constructions The 
idea of twisting [|7|] is to decompose the fields of the Euclidean SYM theory in D spacetime dimen- 
sions in representations not in terms of the original (Euclidean) rotational symmetry SO rot (D), but a 
twisted rotational symmetry, which is the diagonal subgroup of this symmetry and an SOr(D) sub- 
group of the R-symmetry of the theory, that is, SO(D)' = diag(50Lorentz(^) x SOr(D)). It should 
be noted that the R-symmetry group of the theory must be large enough to contain SOr(D) as a 
subgroup. Otherwise, the twisted lattice constructions would not turn out to be successful. 

Two-dimensional £2 = 4 lattice SYM theory: The two-dimensional £2 = 4 SYM theory is 
the simplest example of a gauge theory that permits topological twisting and thus satisfies our 
requirements for a supersymmetric lattice construction. Its R-symmetry possesses an SO(2) sub- 
group corresponding to rotations of its two degenerate Majorana fermions into each other. After 
twisting, we end up with fermionic fields (tj, \p a , Xab) and a set of complexified bosonic fields 
(£/ a , s^a)- The prescription for discretization is somewhat natural. The complexified bosonic 
fields are represented as complexified Wilson gauge links £/ a (x) —> ^a(n), living on the links of a 
lattice, which for the moment can be thought of as hypercubic, with integer-valued basis vectors 
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jUj = (1,0), pL 2 = (0, 1) . They transform in the usual way under U (N) lattice gauge transforma- 
tions %,(n) — > G(n)^,(n)G^(n) . Supersymmetric invariance then implies that y a (n) live on the 
same links and transform identically. The scalar fermion v\ (n) is clearly most naturally associated 
with a site and transforms accordingly T](n) — > G(n)r](n)G t (n) . The two-form field %ab lives on 
the diagonal link. The field orientations are chosen such that gauge invariance is preserved on the 
lattice. 

The continuum derivatives are replaced by difference operators. Motivated by the natural 
technology developed for applying the derivative operators to arbitrary lattice p-forms [§], we need 
just two derivatives given by the expressions: 

^ +) f b (n) = & a (n)f b (n + fi a )-Mn)<& a (n + fi b ) , (2.1) 

^fVa(n) = f a (n)W a (n)-W a (n-iX a )f a (n-ila) ■ (2-2) 

The lattice field strength is then given by the gauged forward difference acting on the link field: 
^afo(n) = f^i + ^%,(n), and is automatically antisymmetric in its indices. Similarly the covariant 
backward difference appearing in @ a ^ a (n) transforms as a 0-form or site field and hence can be 
contracted with the site field tj (n) to yield a gauge invariant expression. 
The final lattice version of the action is 

S = £Tr (j^(n)J^(n) + 1 (^-^(n)) 2 - Xab (n)^ ¥b] (n) - r? (n)^ V«(»)) , (2.3) 

which can be shown to be J2-exact same as its continuum counterpart. 

Four-dimensional £2 = 16 lattice SYM theory: The lattice action of this theory contains a 



jS-exact term of precisely the same form as the two-dimensional theory provided in Eq. ( |L3| ) if 
one extends the indices labeling the fields to run now from one to five. In addition, the appropriate 
twist (called the Marcus twist ^) of jY = 4 YM requires a new i?-closed term, which was not 
possible in the two-dimensional theory, whose lattice version is: 

Wd = - \ £ Tr e mnpqr % qr {n + /l m + 'fi ll + y.p)® ( p\mn (n + ju p ) , (2.4) 

° n 

and can be seen to be supersymmetric since the lattice field strength satisfies an exact Bianchi 
identity [§]: E mnpq ^ ( p ] ^ qr = 0. 

In the case of the jV = 4 theory, the resulting lattice has a nontrivial structure: It is known 
as the A^-lattice. The lattice is constructed from the set of five basis vectors 7 a pointing out from 
the center of a four-dimensional equilateral simplex out to its vertices together with their inverses 
— l? a . Complexified Wilson gauge link variables % a are placed on these links together with their 
cS-superpartners \j/ a . Another 10 fermions are associated with the diagonal links e a +7 b with a>b. 
Finally, the exact scalar supersymmetry implies the existence of a single fermion for every lattice 
site. It is invariant under the exact scalar supersymmetry £2, lattice gauge transformations G(n), 
and a global permutation symmetry S 5 , and can be proven free of fermion doubling problems. The 
J2-exact part of the lattice action is again given by Eq. ( |2T3| ) where the indices a, b now correspond 
to the indices labeling the five basis vectors of the A\ lattice. 

The renormalization of this theory has been recently studied in perturbation theory with some 



remarkable conclusions [13]. 
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Figure 1: The histograms of the phase angles a for £2 = 4, 16 theories with gauge group U(2), A = 1, L = 5. 

Gauge link parameterizations: There exist two distinct parameterizations of the gauge fields 
on the lattice and they both have been used in various simulations of lattice SYM theories. The first 
one follows the standard Wilson prescription where the complexified gauge fields in the continuum 
are mapped to link fields %(n) living on the link between n and in n + ji a through the mapping 
^(n) = e iz/ "( n ), where s^ a {x\) = J\fj si^T 1 and T' = 1,. . . ,Nq are the anti-hermitian generators of 
U (N). The resultant gauge links belong to GL(N,C). We call this realization of the bosonic links 
the exponential or group based parametrization 1 . 

The other parametrization of the bosonic link fields that has been used, particularly in the 
orbifold literature, simply takes the complex gauge links as taking values in the algebra of the U (N) 
group 2 . To obtain the correct lattice kinetic terms needed to target the continuum theory one must 
then expand the fields around a particular point in the moduli space of the theory corresponding 
to giving an expectation value to the trace mode of the imaginary part of the link field - a field 
which can be identified as a U(l) scalar field in the untwisted theory. The expansion is: ^,(n) = 
J^+«fz&(n). Usually the use of such an algebra based or non compact parametrization would signal 
a breaking of lattice gauge invariance. It is only possible here because the bosonic fields take values 
in a complexified U (N) theory - so that the unit matrix J? appearing in the expansion above can 
be interpreted as the expectation value of a dynamical field. We will refer to this parametrization 
as the linear or algebra based parametrization. 

Both parameterizations of the gauge links are equivalent at leading order in the lattice spacing, 
yield the same lattice action and can be considered as providing equally valid representations of the 
lattice theory at the classical level. In this work we have concentrated on the linear parametrization 
principally because it is naturally associated with a manifestly supersymmetric measure in the path 
integral - the flat measure. 

The correct continuum limit: The requirement that the theory target the correct continuum 
theory requires that the fluctuations of all dimensionless lattice fields should vanish in the contin- 
uum limit. In addition the linear parameterization only yields the correct naive continuum limit if 
the trace mode of the scalars develops a vacuum expectation value so that appropriate kinetic terms 
are generated in the classical action. Since no classical scalar potential is present in the lattice 
theory it is crucial to add by hand a suitable gauge invariant potential to ensure this feature. 

Specifically we add a potential term of the form Sm = A* 2 Lx {jj^r^J (x) c 2^ a (x)) — l) 2 to the 
lattice action Here }X is a tunable mass parameter, which can be used to control the expectation 

1 Notice that our lattice gauge fields are dimensionless and hence contain an implicit factor of the lattice spacing a. 

2 In fact, a non-compact parametrization of the gauge-fields is also recently used to restore the BRST symmetry on 
the lattice in Ref. Jl(i|], i.e., to evade the so-called Neuberger 0/0 problem Jll| ] (see also Refs. jlo| and [ |l2| for the 
recent progress). 
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Figure 2: The Pfafian phase < | sin a \ > for =2 = 4,16 theories with gauge group U (2) and ji =0.1,1,10. 

values and fluctuations of the lattice fields. Notice that such a potential obviously breaks super- 
symmetry - however all supersymmetry breaking counter terms induced via quantum effects will 
possess couplings that vanish as pL — > and so can be removed by sending jU — > at the end of the 
calculation. 

In our simulations, we have rescaled all lattice fields by powers of the lattice spacing to 
make them dimensionless. This leads to an overall dimensionless coupling parameter of the form 
N/(2Xa 2 ), where a = p/T is the lattice spacing, j3 is the physical extent of the lattice in the Eu- 
clidean time direction and T is the number of lattice sites in the time-direction. Thus, the lattice 
coupling is K = , with t = Aj3 2 , for the symmetric two-dimensional lattice where the spatial 
length L = T . Note that t is the dimensionless physical 't Hooft coupling in units of the area. In our 
simulations, the continuum limit can be approached by fixing t and N and increasing the number 
of lattice points L — > oo. We have taken three different values for this coupling t = 0.5, 1 .0, 2.0 and 
lattice sizes ranging from L = 2, ■ • • ,12. We fix the value of j8 = 1 so that t = X. Theories with 
U (N) gauge group with N = 2 have been examined. 

The simulations are performed using anti-periodic (thermal) boundary conditions for the fermions 
(This forbids exact zero modes that are otherwise present in the fermionic sector). An RHMC al- 
gorithm was used for the simulations as described in [[Tall- 
in Fig. (|I]), we show histograms of the phase angles a for different L and /I = 0.1, 1, 10 for 
both <£ = 4 and £2 = 16 theories (for A = 1 and 0.5, respectively) with gauge group U(2). In 
Fig. (|J) we show the results for the absolute value of the (sine of) the Pfaffian phase 3 . Three values 
of /j. are shown corresponding to ju = 0.1, 1, 10. While modest phase fluctuations are seen for small 
lattices for the case of jj. = 0.1 we see that they disappear as the continuum limit is taken. 



3. Conclusions 

We have performed numerical simulations of the four and sixteen supercharge lattice SYM 
theories in two dimensions to investigate the occurrence of a sign problem in these theories. In 
contrast to the usual situation in lattice gauge theory, we utilize a non compact parameterization 
of the gauge fields in which the lattice fields are expanded on the algebra of the group. We have 
examined both supersymmetric theories for several values of the dimensionless 't Hooft coupling 
A/3 2 and for the gauge group U (2). We take a careful continuum limit by simulating the theories 

3 Obviously, the absolute value of sine of the angle does not distinguish between, e.g., and n. However, focusing 
on a single measure of the magnitude of the phase fluctuations does allow for a clean extrapolation of to the continuum 
limit and, as our previous histograms show, suffers from no ambuiguity. 
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over a range of lattice size L = 2 — 12. In both cases we observe that the average Pfaffian phase 
is small and goes to zero for a fixed gauge invariant potential as the continuum limit is taken. 
Indeed, in practice it is sufficiently small even on coarse lattices that there is no need to use a 
reweighting procedure to compute expectation values of observables. The detailed analysis and 
results (including for the U (4) theories) are provided in [|i"5|]. 
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